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^ ■ ABSTRACT 

O ! I present sensitive high resolution VLA B, C, and D array observations of the almost 

edge-on Scd galaxy NGC 4244 in the 21-cm spectral line of neutral atomic hydrogen. 
The gas layer of NGC 4244 is rather symmetric in all respects, i.e. the surface density 
distribution, flaring and warping. This symmetry allows for a reliable determination of 
the rotation curve, despite the fact that the galaxy is close to edge-on. The rotation 
curve rises slowly in the inner 6 kpc, is roughly constant at 100 km s~^ out to 10 kpc, 
^ . and decreases in Keplerian fashion by 15% at the last measured point at 14 kpc. The 

rotation curve constrains the stellar mass-to-light ratio to lie between 50 and 100% of 
the "maximum- disk" value. 
^ ' A new technique is presented to determine simultaneously the inclination and the 

. thickness of the gas layer from high resolution H I observations. This procedure uses 

the apparent widths at many azimuths (many channels) and can be used at inclinations 
as low as 60° . Kinematic information is used to separate flaring from warping. The 
inclination of the unwarped disk is about 84?5, while the small warp coincides with a 
decreasing inclination (to 82?5 ± 1° ). The data indicates that at large radii the disk 
warps back to the plane defined by the inner disk. The measured gaseous velocity 
dispersion is roughly constant within the optical disk (8.5 ± 1 km s~^) and increases 
slightly beyond. 

On both sides of the galaxy the thickness of the gas layer increases gradually from 
~400 pc at 5 kpc to ~1.5 kpc at the last measured point (at 13 kpc). The strong 
gradients in the inferred thickness which bracket the spiral arms probably result from 
streaming motions associated with the arms and are not intrinsic to the galaxy. In an 
accompanying paper (AJ, Aug. 1996) I use the measurements presented in this paper 
to infer that the dark matter halo of NGC 4244 is highly flattened. 
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1. INTRODUCTION 

In a previous paper, Paper I (Oiling 1995), 
it was shown that measurements of the radial 
variation of the thickness of gas layers (flaring) 
can be used to constrain the shape of the dark 
matter (DM) halo. This is accomplished by com- 
paring the measured flaring with that expected 
from a self-gravitating gaseous disk in a poten- 
tial generated by the stellar disk and (flattened) 
DM halo. Here, high resolution H I observation 
of the nearby, almost edge-on, Scd galaxy NGC 
4244 are presented. Using these observations I 
determine the rotation curve, the flaring, the ve- 
locity dispersion, and the radial surface density 
distribution of the gas. In the accompanying pa- 
per (Paper III, AJ, 1996, ...) I will apply the 
method presented in Paper I to determine the 
flattening of the dark matter halo of NGC 4244 
using the essential parameters as derived in this 
paper. 

In a study similar to the present one, Rupen 
(1991a) analyzed high resolution H I observations 
of two edge-on galaxies, NCC 4565 and NCC 
891. Due to the increased sensitivity of the 
VLA, sensitive high resolution H I observations, 
as pioneered by Rupen (1991b), can now be rou- 
tinely obtained and processed. The high reso- 
lution H I data presented in this paper are sen- 
sitive enough to determine the gaseous velocity 
dispersion and probe the region beyond the op- 
tical disk where the thickness of the gas layer is 
most sensitive to the shape of the DM halo. The 
new technique to determine the thickness of the 
gas layer (presented in Sec. |5]) complements an 
existing procedure specific to the case of strong 
spiral structure in the H I (Braun 1991), and is 
more general than the method employed by Irwin 
& Seaquist (1991). 

NGC 4244 is an intermediate mass system 
with a rotation speed of about 100 km/s. Its 
proximity (3.6 Mpc) results in the high linear 
spatial resolution (17.5 pc/arcsec) which is nec- 
essary to successfully determine both the verti- 



cal structure of the gaseous disk and the rota- 
tion curve in the inner parts of the galaxy. The 
large inclination of NGC 4244 with respect to the 
line of sight, about 84.5 degrees, in combination 
with the symmetric surface density distribution, 
allows for an accurate determination of the rota- 
tion curve. NGC 4244 is one of the few galaxies 
to show, on both sides, the characteristic drop in 
rotation velocity beyond the optical disk. 

The outline of this paper is as follows: in 
Sec. ^ I summarize current observational knowl- 
edge of NGC 4244. The H I observations are 
presented in Sec. ^ and the parameters derived 
therefrom in Sec. 0. A new technique to derive 
the thickness and inclination of the gas layer is 



outlined in Sec. 5.1., and is applied in Sec. 5.2 



The results are discussed in Sec. ^ 

2. THE GALAXY NGC 4244 

NGC 4244 is a small almost edge-on ScdIII 
galaxy in which individual stars and H II regions 
can be resolved (Sandage & Bedke 1994). Sur- 
face brightness photometry (van der Kruit & 
Searle 1981a) reveals that the radial light dis- 
tribution is close to exponential till about five 
optical scale lengths, where it cuts off sharply 
(see also Fig. Q). The vertical scale height of 
the light is constant with radius, while no thick 
disk is present. The very faint bulge (/iK('^e) = 
23.05 mag/arcsec^) with an effective radius of 
140" (Bergstrom et al. 1992) is dynamically in- 
significant. 

There are several indications that the cur- 
rent star formation rate is very low in NGC 
4244. Only very low radio continuum emission 
(Condon 1987) and little Ha emission (Walter- 
bos 1994, private communications) have been de- 
tected, while upper limits have been set on X- 
ray emission (Bregman & Glassgold 1982) and 
emission from highly ionized halo gas (Dehar- 
veng et al. 1986). Very faint radio continuum 
emission extends along the plane of the galaxy 
from 2.'5 north-east (NE) to 2' south-west (SW) 
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of the center. The 60/xm emission is enhanced 
in this region as weh (Rice 1993). In contrast to 
most galaxies (Bicay & Helou 1990), NGC 4244's 

20 cm radio continuum emission is not more ex- 
tended than the 60/xm emission. The weak CO 
emission (Sage 1993), coincident with the radio 
continuum emission, amounts to 2.2 x 10^ Mq of 
H2 (50% of the H I mass in this region). The H2 
radial surface density distribution is calculated 
with the Abel inversion technique (Sec. |4.1.| ) and 
is consistent with a Gaussian ring centered at 1.3 
kpc, with a FWHM of 1.6 kpc and an amplitude 
of 1.6 M0PC-2. 

Throughout this paper the distance to NGC 
4244 is taken to be 3.6 Mpc, consistent with 
the IR-Tully-Fisher distance (3.7 Mpc, Aaron- 
son et al. 1986) and the "hydrogen counterpart" 
of the Tully-Fisher distance (3.6 Mpc) as defined 
by Broeils (1992, Chapter 3). NGC 4244 is the 
second brightest member of a group of dwarf 
galaxies: the CVn I (de Vaucouleurs 1975) or B4 
(Kraan-Korteweg & Tammann 1979) group with 

21 catalogued members. This group has a ra- 
dius of 1 Mpc, an average density of 4.5 Mpc~^, 
and a total luminosity of 8.9 x 10^ Lq^b- The 
nearest catalogued neighbor is the SmIV galaxy 
NGC 4190 at a projected distance of 81 kpc. This 
group is rather elongated (axial ratio ~ 3), with 
the position angle of the long axis at -18° , 66° 
from the plane of NGC 4244. The observation 
presented in this paper (Sec. |3.1.|) fail to detect 
gas rich companions more massive than lO"^ Mq 
within 45 kpc. 

Figure |l| shows the optical (grey scale) image 
from IllaJ photographic observations (van der 
Kruit & Searle 1981a) and a total hydrogen map. 
The orientation of the weak dustlane suggests 
that the eastern side of the galaxy is the far side 
and thus least affected by absorption. The intrin- 
sic radial surface brightness distribution, L{R), 
is derived by applying the Abel-inversion tech- 
nique (Sec. 4.1.) to four different strips parallel 



dicular to the minor axis (top profile in Fig. ^). 
The second strip contains all light within 420 pc 
(twice the vertical exponential scale height, Zg) 
from the midplane. The third and forth strips 
are determined from those parts which are fur- 
ther than 420 pc from the midplane (east and 
west respectively). As the region which is af- 
fected by dust extinction lies generally close to 
the plane {z <^ Vs'Ze ;Peletier et al. 1995), the 
third and forth strips are expected to be free 
of extinction effects. All four profiles indicate 
a radially exponential light distribution with the 
same radial scale length (/ir) of 2.0 kpc, slightly 
larger than the value derived by van der Kruit &; 
Searle (1981a)Q. Thus, the overall light distribu- 
tion is probably not much affected by extinction 
(the adopted fit is indicated by the filled circles). 
As noticed by van der Kruit & Searle (1981a), the 
stellar disk is truncated at somewhat different ra- 
dial distances on the two sides of the galaxy. The 
IllaJ-magnitudes are converted to B-magnitudes 
assuming a (B-V)=0.7 as appropriate for the old 
disk (Bottema 1989). Some parameters of this 
galaxy are listed in Table 

3. H I OBSERVATIONS 

NGC 4244 was observed with the VLA| 
(Napier et al. 1983), in the 10 km (B), 3 km (C), 
and 1 km (D) arrays for a total of about 14 hours. 
The data from 1989 and 1990 were taken as part 
of a project aimed at doing a multi-wavelength 
high resolution study of the interstellar medium 
in nearby galaxies (Braun 1995). The observing 
parameters are summarized in Table ^. 

Standard calibration procedures were followed. 
After removal of some data affected by inter- 



to the major axis, after removal of foreground 
stars. The first strip contains all light perpen- 



^AU the linear distances derived by van der Kruit & Searle 
(1981a), who use a distance of 5 Mpc (Sandage & Tam- 
mann, 1975), are scaled to the adopted distance for NGC 
4244, 3.6 Mpc. 

^The VLA of the National Radio Astronomy Observatory 
is a facility of the National Science Foundation operated 
under cooperative agreement by Associated Universities, 
Inc. 
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ference^ the data were combined in the UV plane. 
The primary flux calibrator was 3C 286, with an 
assumed flux at 1418 MHz of 14.86 Jy. A heav- 
ily tapered continuum image was made using 5 
line-free channels on both sides of the band. The 
rms noise in this image is 0.94 mJy beam~^ at 
a resolution of 80" x 77" . Several background 
sources but no extended radio continuum emis- 
sion was detected from NGC 4244 itself down 
to a 3 sigma limit of 2.8 mJy beam~^, consis- 
tent with a peak flux of 1.4 mJy as determined 
by Condon (1987) and Hummel et al. 's (1984) 
non-detection at 610 MHz. A strong background 
source about 7' NE of NGC 4244 is probably the 
reason for the much stronger detection reported 
by Schlickeiser et al. (1984). 

For the two high resolution data sets the sen- 
sitivities listed in Table || correspond to roughly 
3 X lO'^ Mq beam~^ channel"^ at the field cen- 
ter, while the two low-resolution sets are about 
three times more sensitive as a result of the in- 
terim upgrade of the VLA receivers. To bring 
out emission at different scales, images of the line 
emission were made with different tapering. 

3.1. Results 

Assuming low optical depth the total amount 
of neutral hydrogen is arrived at by adding the 
flux containing areas in the channel maps. The 
validity of this assumption is checked in Ap- 
pendix The flux containing areas are defined 
to be those regions which have a flux larger than 
^ 1.2a in a smoothed version of the spectral 
line cube (smoothed to a four and three times 
lower resolution, spatially and in velocity res- 
pectively). The high and low-resolution total 
H I distributions are presented in Fig. The 
images have been rotated 48° clockwise, so that 
the upper part is the north-eastern side. Both 
the high and low resolution data show a warp 
starting at the edge of the optical disk (about 9 

■^In particular the 1992 D-array daytime observations re- 
quired extensive UV-editing due to solar interference. 



arcmin from the center). The H I distribution of 
NGC 4244 is rather symmetric apart from the 
"popped blister" located just below the major 
axis about four arcmin to the north-east. A re- 
cent star formation event might be responsible 
for this feature as it coincides with peaks in the 
radio continuum and IRAS emission. 

The low resolution, primary beam corrected, 
total intensity image yields a total H I mass of 
(1.29 ± 0.05) X 10^ Mq. The north-eastern part 
contains 51% of the total H I mass, the south- 
western part amounts to 49%. The hydrogen 
mass determined from the global proflle (the flux 
as a function of velocity. Fig. |), (1.4±0.1) x 10^ 
Mq, is consistent with the determination from 
the total H I map. The H I mass listed in Table 
^ is the average of the total H I and global profile 
values, where the error equals half the difference. 

From the global proflle a systemic velocity 
(the mean of the 10%, 20% and 50% points) 
of 244.4 lb 0.3 km is determined, in agree- 
ment with the value determined by Huchtmeier 
& Seiradakis (1985). As a result of the asym- 
metry in the H I distribution, the receding part 
has a higher peak value than the approaching 
part. The position at the major axis where the 
recession velocity equals the systemic velocity is 
offset by 7"4ib0"6 towards the south-west, which 
agrees well with the central position of the H I as 
determined in Sec. |4.1. . 

Within 45 kpc from the pointing center, no 
companions were detected down to a flve sigma 
level of 10^ Mq channel"^ at the center (or 1.6 x 
10^ Mq channel"^ at 45 kpc). 

Fig. 1^ shows a composite of the high resolu- 
tion (28"3xl0"l or 490x177 pc along the major 
and minor axes respectively) and intermediate 
resolution (38"x38") channel maps. The maps 
are combined such that the channels with the 
same speed relative to the systemic velocity are 
present in one frame (after having blanked those 
areas without emission). The receding part of 
NGC 4244 is the south-western part (negative 
X- values), the north-eastern part is approaching. 
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The low velocity channels, which would overlap 
if presented in this manner, are presented indi- 
vidually. One local irregularity is clearly visible: 
the +4' "hole" found in the total H I map can be 
seen in the frames labeled 17,18, and 19. 

Each map shows a characteristic "filled V" 
shape. Two different mechanisms can cause such 
shapes. First of all it reflects the flaring of 
the gaseous disk. The other contributing fac- 
tor is kinematic in nature: the "butterfly" pat- 
tern (seen in channel maps of less inclined sys- 
tems, e.g., Bosma 1981) seen almost edge-on. In 



Sec. 5.1. I present a new technique to disentan- 
gle these two effects. Classically the flaring is 
measured in the edge channels, i.e. at velocities 
of {Vsys ± Vrot,max) (Sancisi & Allen 1979; van 
der Kruit 1981; Rupen 1991a). Because NGC 
4244 is not exactly edge-on and because the ve- 
locity edge depends on radius (due to the falling 
rotation curve), the apparent widths can not be 
easily interpreted in terms of an intrinsic verti- 
cal scale height. Furthermore, systematic errors, 
e.g. due to non-circular motions are largest in the 
edge channels (Appendix^]). The channels ±18 
are the edge channels between about 5' and 9' . 
For this radial range, the receding side exhibits a 
more-or-less regular increase of apparent width, 
but on the approaching side this effect is not ev- 
ident. This difference can easily arise from local 
effects such as arm inter-arm surface density dif- 
ferences, holes, bubbles, and streaming motions. 
The new thickness-and-inclination determination 



technique outlined in Sec. 5.1. minimizes the ef- 
fects of such irregularities. 

A rather regular warp, starting just before the 
optical edge of the galaxy, is seen in the channel 
maps. The channels ±16, ±17, and ±18 show the 
strongest warping at large radii. The position 
angle of the warp is determined by tracing the 
major axis of emission as a function of major 
axis distance. Within the optical disk, the posi- 
tion angle of the centroid is constant at 48° east 
of north. At a radius of 8' , the position an- 
gle begins to change linearly to reach 43° at 13' . 



Because no clear butterfly pattern is seen in the 
warped region, the warp is likely to be the result 
of a change in position angle rather than in incli- 
nation. This is confirmed by the derived radial 
variation of the inclination (Fig. ^). 

4. INFERRED PROPERTIES 

In this section several properties of the gas^ 
layer will be determined: the radial surface den- 
sity distribution (Sgas(-R)), the rotation curve 
{Vrot{R)), and the radial variation of the gaseous 
velocity dispersion ((Tgas(-R))- All these quanti- 
ties are required to make self-consistent predic- 
tions of the radial variation of the thickness of 
the gas layer ( Paper I, Paper III). Furthermore, 
with the aim to gauge the reliability of the de- 
termination of these properties, these quantities 
are used to compare the VLA observations with 
simulated observations of NGC 4244 (Appendix 



Integrating the total hydrogen along the mi- 
nor axis results in a major axis H I profile (A(x)). 
The inner 4' of this profile are symmetric about a 
point located 6".9 ± 0''6 to the south-west (along 
the major axis) of the catalogued position. Sim- 
ilarly, the location of the midplane is offset by 
2'f2 ± 0".5 to the south-east. In the rest of this 
paper I assume that NGC 4244's center is lo- 
cated at the centroid of the H I -distribution: a 
= 12^ 15™ 9^7, 6 = 38° 4' 59"A, which position is 
listed in Table ^. This position agrees with the 
kinematic center determined in Sec. |3.1.| . 

4.1. The H I Surface Density Distribution 

I describe the line of sight deconvolution nec- 
essary to determine the radial surface density 
distribution after introducing the coordinate sys- 
tems used. First there is the coordinate system 
centered on the center of the galaxy: (x, y, z) in 



*Unless otherwise noted, the term "gas" refers to H I in an 
observational context, and to the sum of H I , H2 and He 
in a dynamical context. 
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Cartesian coordinates and {R, 6, z) in polar co- 
ordinates (with = being the intersection of 
the plane of the sky with the plane of the ga- 
laxy, and with z the vertical coordinate). In the 
observers coordinate frame, x and d are coinci- 
dent with the major and minor axes of the inner, 
unwarped, part of the disk. The projection of 
the line of sight onto the (unwarped part of the) 
plane of the galaxy is parallel to the y-axis. 

I determine an approximation to the H I profile 
(the total H I distribution summed perpendicu- 
lar to the major axis) by summing the emis- 
sion over the d-axis. Because of the warp the 
x-coordinate does not truly trace the major axis 
(x'). However, since x' = ^/x"^ + (P does not dif- 
fer significantly for the small warp of NGC 4244 
{d < 1 kpc for X < 13 kpc), the derived density 
distribution is not affected. The integration^ of 
the volume density along the line of sight and 
perpendicular to the major axis, which yields the 
H I profile, can be written as an integration of the 
true surface density along the y-axis: 



Hx') 



2 X 



dy Sgas 

R Sgas(^) 
- x2 



x"^ + y2 



dR (Mo/pc). (1) 



Inversion of this Abel integral [e.g., Binney & 
Tremaine, Eq. (lB-59b)] yields the radial surface 
density distribution 



Sgas(-R) 



1 



dK{x) 



In (Me/Pc=).(2) 

This procedure can be applied at any wavelength 
for any inclination provided that the optical 
depth is low, the galaxy is not warped and if it is 
axisymmetric. Note that similar techniques were 
used to determine the surface brightness distribu- 
tion of Virgo spirals from "pencil-beam" H I data 
(Warmels 1988) and the mass-to-light ratio of el- 
liptical galaxies (i.e., Schwarzschild 1954; Binney 



''Assuming low optical depth, see Appendix 



&; Mamon 1982). Adding realistic noise to a sim- 
ulated spectral line cube (Appendix B.) and ap- 
plying Eq. (|2|) yields the input Sgas(^). Thus, 
the derivation of the surface density distribution 
is not sensitive to noise. However, if the signal to 
noise ratio is poor, it may be preferable to use the 
iterative technique proposed by Warmels (1988). 

5^gas(-R) is derived from a combination of the 
high and low resolution A(x) curves and is pre- 
sented in Fig. ^ Negative surface densities 
at small radii occur unless the two sides of the 
A(x) profile are averaged over the inner ~ 1.5 
kpc. Thus this method does not recover the true 
Sgas(-R, 0) but it can find an azimuthally sym- 
metric density distribution consistent with the 
observations: it serves as a reasonable approxi- 
mation to the true distribution. Integrating the 
E(i?) profiles out to 15.8 kpc (where the surface 
density equals 0.02 Mopc'^), I find 1.4 x 10^ M©: 
the Abel inversion technique conserves flux. 

Several strong (up to a factor 3) density en- 
hancements (~1 kpc wide), possibly spiral arms, 
are present. In the inner 4 kpc the south-western 
and north-eastern density distributions are very 
similar. Beyond 10 kpc the surface density de- 
creases in an exponential manner with a scale 
length of roughly 1.36 kpc. No sharp cutoff in 
the H I column density is seen above the 3cr de- 
tection limit (5 X 10^^ cm~^). We see that the 
surface density distribution is rather symmetric 
(see also Sec. |5.2.| ), except in the 6-8 kpc re- 
gion. It might be that this surface density dif- 
ference causes differences in the observed flaring 
(Fig. 13) as a result of the different self- gravity 
at these points. 

4.2. The Rotation Curve and the Veloc- 
ity Dispersion of the Gas 

Figure |^ shows the intensity as a function 
of velocity and position along the major axis 
(summed over the minor axis). The rotation 
curve and velocity dispersion of NGC 4244 are 
derived from this plot. The maximum relative ve- 
locity outline corresponds roughly to the rotation 
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speed minus twice the velocity dispersion in the 
tangential direction. Approximations to the ro- 
tation speed and velocity dispersion were deter- 
mined by fitting a single Gaussian to the extreme 
velocity envelope. These approximations have 
been corrected for instrumental, line-of-sight in- 
tegration, and beam smearing effects, as detailed 



in Appendix A. 

Figure ^ shows NGC 4244's smearing cor- 
rected intrinsic tangential dispersion as derived 
from the high, intermediate and low resolution 
data sets, which sample different ranges in galac- 
tocentric radius. For each set of measurements, 
both sides yield dispersions generally equal to 
within the errors. The errors on the values ex- 
tracted from the lO'f data set are large because 
the signal to noise ratio is relatively low. The av- 
erage dispersion curve is the weighted average of 
the measurements obtained at the two sides (38" 
data set for R £ (9,13] kpc, and 55" beyond). 
For radii smaller than 9 kpc I also average the 
measurements from the Hanning-smoothed 10"l 
and 38" data sets since these data arise from es- 
sentially independent measurements (B & D ar- 
ray, respectively). 

In Fig. ^ I present the line-of-sight and beam- 
smearing corrected rotation curves for both sides 
of the galaxy. Apart from the region between 
2.5 and 7.5 kpc, the rotation curves are equal to 
within the errors. The drop in the north-eastern 
rotation curve seen between 2.5 and 5.5 kpc is 
clearly a result of a local phenomenon (described 
above): there the adopted rotation curve follows 
the south-western curve (10" xlO" ) till 4.5 kpc. 
The adopted rotation curve between 4.5 and 6 
kpc is a smooth junction of the adjacent regions: 
here the rotation speeds as determined from the 
XV plots are probably influenced by the strong 
H I peak at 5.5 kpc (see Appendix ^ for a dis- 
cussion). The 6 to 7.5 kpc region, on the south- 
western side, corresponds to the region where the 
B-band surface brightness is depressed by a fac- 
tor of about 10. Therefore, I choose the intrinsic 
rotation curve to follow the north-eastern curve 



in this range of galactocentric radii (38" x38" ). 
Taking these modifications into account leads to 
the adopted intrinsic rotation curve (indicated by 
the filled circles in Figure ^), which is minimally 
affected by non-circular motions. 

Because of the large inclination, small changes 
thereof are not important for the derivation of 
the rotation curve. The maximum amplitude of 
the rotation velocity is about 100 km s~^ and 
starts dropping significantly beyond -Rmaxj the 
edge of the optical disk. At a distance of 14 kpc 
the rotation speed has dropped to 85 kms~^. 
Huchtmeier (1975) finds, using a 9' beam, that 
the H I extends out to 18' and determines a ro- 
tation speed of 85 lb 5 km s^^ at that distance. 
The measurements presented in this paper do not 
confirm this detection. The symmetry between 
the two sides indicates that the gas is likely to 
be on circular orbits, except for the regions men- 
tioned above. The symmetric drop in rotation 
curve on both sides of the disk is the strongest 
argument that we are actually measuring the cir- 
cular velocity and that this drop is not caused by, 
for example, a lack of gas on the major axis: an 
explanation which has been invoked to explain 
asymmetric drops seen in other edge-on galaxies 
(Sancisi & Allen 1979). 

5. THE DETERMINATION OF THE 
THICKNESS OF THE HYDRO- 
GEN LAYER 

In the past, the radial increase of the thickness 
of the gas layer of external galaxies has been de- 
termined from apparent width measurements of 
the edge channels only (Sancisi & Allen 1979; van 
der Kruit 1981; Rupen 1991a). The main advan- 
tage of this procedure is that the edge channels 
select a region (the "emitting region" ) around the 
major axis which is relatively narrow, so that its 
width, as projected onto the sky, does not con- 
tribute much to the observed apparent thickness. 
This procedure has two disadvantages however: 
1) the magnitude of the projection effect depends 
strongly upon inclination, e.g., at i=85° it varies 
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by a factor of two for an inclination difference of 
2° , and 2) only two points per galactocentric ra- 
dius, one on the approaching side and one on the 
receding side, axe measured. 

Irwin k, Seaquist (1991) modeled the entire 
spectral line cube obtained for NGC 3079 to ob- 
tain the rotation curve, surface density distribu- 
tion and gas layer width (assumed to be constant 
with radius). Braun (1991) devised a technique 
to determine the inclination and the thickness of 
the gas layer for M31, where the H I emission is 
found primarily in the spiral arms. The tech- 
nique presented here works best if the H I is dis- 
tributed smoothly. 

Below I describe a new technique to simul- 
taneously determine the radial variation of the 
inclination and the thickness of the gas layer of 
external galaxies. As the essence of this method 
is to utilize data from many azimuths (channels) , 
both disadvantages of the edge-channel proce- 
dure are circumvented. This is accomplished by 
modeling the contribution of the emitting region 
to the apparent width. Two aspects are crucial: 
1) the emitting region can be well modelled if the 
rotation curve is known, for example for a flat ro- 
tation curve a channel of a several km s~^ wide 
is sensitive to gas in a wedge-like region in the 
plane of the galaxy (Sancisi & Allen 1979, their 
Fig. 6), where the width of this emitting region 
follows directly from the rotation curve, the ob- 
serving velocity and velocity resolution, and 2) 
for a given galactocentric radius the planar width 
of the emitting region varies systematically with 
azimuth (i.e., observing frequency, or velocity). 
Thus the projected width of the emitting region 
also varies systematically with azimuth and ra- 
dius in a manner which can be described exactly 
for a given observed rotation curve and channel 
width. These projection effects can be clearly 
seen in Fig. 8 of Sancisi & Allen, where the ap- 
parent width of the gas layer is larger in those 
channels whose velocities are closer to the sys- 
temic velocity. The apparent width of the gas 
layer is then the sum of the projected vertical 



thickness of the gas layer and the projected pla- 
nar width of the emitting region. Only two un- 
knowns (the vertical thickness and the inclina- 
tion) are involved. Since the apparent widths 
can be measured at many azimuths, and the pla- 
nar width of the emitting region can be calcu- 
lated as a function of galactocentric radius and 
azimuth, the inclination and vertical thickness of 
the gas layer can be simultaneously determined 
from measurements of the apparent thickness in 
many channels. 

5.1. Outline of the Procedure 

The procedure outlined below works for a 

more realistic situation where the rotation curve 
is not flat and where the gaseous velocity disper- 
sion softens the edges of the wedges. A more 
detailed description is given elsewhere (Oiling 
1995b), and is applied in the next subsection. 

In the case where a galaxy is not edge-on, 
the azimuthal dependence of the flaring might 
be determined from measurements of the appar- 
ent thickness (in a given velocity channel) along 
an ellipse with axial ratio 1/ cosi. For a close to 
edge-on galaxy this procedure is difficult, while 
at z = 90° it fails altogether. However, at any 
inclination, the galactocentric radius {R) can be 
determined from the rotation curve (T4ot(-R))) the 
major axis position {x) and the observing veloc- 
ity (Khan): 



R 



V,ot{R) smi{R) 



cos 9 



Vr 



chan 



sys 



X 



Khan — Vs 



(3) 



sys 



with i(R) the inclination, and Vgyg the galaxy's 
systemic velocity. Here the assumption is made 
that the galaxy can be described adequately by a 
tilted ring model (Rogstad et al. 1974; Begcman 
1989) which implies that the warp angle {i{R)) 
does not depend on azimuth. It is important 
to realize that the determination of R depends 
on the projected rotation speed (i.e. the product 
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of VrotiR) and s'mi{R)) which is, at any inchna- 
tion, much better determined than the individual 
quantities (Begeman 1989): the kinematic deter- 
mination of R is robust]^. When the apparent 
widths in a given channel are determined, Eq. 
(^) is used to assign a galactocentric radius to 
the measurements obtained at major axis posi- 
tion X. For a given channel and major axis po- 
sition, the planar width of the emitting region 
can be determined as well, so that the apparent 
width measurements can be corrected for inclina- 
tion and projection effects to yield the intrinsic 
thickness of the gas layer. Below I describe in de- 
tail how the width of the planar emitting region 
as well as the inclination and thickness of the gas 
layer can be determined. 

Rather than using the wedge-approximation 
discussed in previous sections, I include the ef- 
fects of the velocity dispersion of the gas, the ex- 
act shape of the rotation curve, and the gaseous 
surface density distribution to determine the pla- 
nar width of the emitting region. Given these 
parameters it is easy to calculate, for all points 
in the galaxy what fraction of the total 

column density can be detected when ob- 

serving the galaxy at inclination i{R) and veloc- 
ityQ Khan (Oiling 1995b, Chapt. 4). These chan- 
nel surface density distributions (CSDDs) form a 
three dimensional data set y, Khan))) where 
each plane represents the face-on surface density 
distribution which has (given the galaxian prop- 
erties and the observing geometry) the right ra- 
dial velocity to be observed in the channel cen- 
tered on Khan- If observed face-on for example, 
the emission observed in each channel is a frac- 
tion (which depends on the gaseous velocity dis- 
persion and the channel response function, but 
is constant across the galaxy) of the surface den- 
sity. Some CSDD maps are presented in Fig. p^ . 
Alternatively, one could think of a CSDD as a 



^The effects of non-circular motions are discussed in Ap- 
pendix 

^The Channel Response Function appropriate for these 
VLA observations is used (see Appendix A]). 



deprojected (to face-on view) channel map of an 
infinitely thin gas layer. A rough idea about the 
shape of the CSDDs can be formed by depro- 
jecting (to face-on view) a velocity field which is 
contoured at Mj4X(AKhan) cTgas) km s~^ inter- 
vals and selecting the area between two contours 
(with AKhan the width of the channels). But 
note that a CSDD is neither a deprojected chan- 
nel map nor a deprojected velocity field. Along 
a line perpendicular to the major axis, the shape 
of the emitting region is no longer discontinuous, 
but can be parameterized, to reasonable accu- 
racy, by two Gaussians: 



^c{y;X, Khan) 



oc 



iiy-Yc)/Wcf 



+ 



e-hi(y+Yc)/Wcf 



(4) 



one for the near side and one for the far side. 
Non-linear least squares fits to the profiles ex- 



tracted from Fig. 10 yield typical errors of 0.5% 
- 1.5% and 0.2% - 0.6% for Wc and Yq, respect- 
ively. As the CSDDs are different for each chan- 
nel, so are (for a given galactocentric radius) Yq 
and Wc, which depend mainly on the rotation 
curve and the gaseous velocity dispersion, res- 
pectively. Projecting the CSDD onto the sky, 
the apparent planar width is given by: 

cos i FWHMp{x; Vchi^n) ~ 2.35 cos i x 

^JW^{X■, Khan) + [Yc{x; Khan)]' , (5) 

where the Yc-term is only included when the pro- 
jections of the near and far side of the CSDD 
are unresolved by the beam, i.e. in nearly edge- 
on viewing geometry. Since the vertical density 
distribution of the gas is approximately Gaus- 
sian (Paper I) , and the planar distribution of the 
CSDD resembles a double Gaussian (Eq. [Q|), 
the line-of-sight integration through this density 
distribution is approximately double Gaussian as 
well. In nearly edge-on geometry the two Gaus- 
sians merge (the butterfly 'closes its wings'). The 
observable width, convolved with the synthesized 
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beam {FWHM\^), is given by: 

FWRM^MVch^n) =FWHMl + 
[FWHM^{R)sini{R)f + 
[FW H Mp{R; Vch.n) cos i{R)f , 



(6) 



Using this equation, a fit to the measured appar- 
ent widths in many channels yields the intrin- 
sic thickness and the inclination of the gas layer. 
The required values for FWHMp{R;Vchan) are 
determined from the CSDD maps. In each chan- 
nel map (and CSDD map) a major axis posi- 
tion is chosen such that galactocentric radius R 
is sampled [via Eq. I solve Eq. (^) us- 

ing a non-linear least squares technique^ which 
also yields the errors on i and FWHM^- For the 
case of NGC 4244, typically 12 determinations of 
FWHMohsiR) are used. The square of the differ- 
ence between the observed and model widths is 
minimized by varying the inclination and gaseous 
width, where each data point is weighed by the 
square of the inverse of its error. The error 
on each data point is calculated by adding (in 
quadrature) the error in the right-hand side of 
Eq. (P), as calculated from the errors in i and 
FWHAIp, to the observational error. Since the 
contribution from the right-hand side depends on 
the exact value of the inclination, several itera- 
tions are required. 

Previous workers (Sancisi & Allen 1979; van 
der Kruit 1981; Rupen 1991a) corrected their 
edge channel flaring determinations for some of 
these projection effects, but made no effort to 
model the emission in lower velocity channels 
in detail. However, as shown above, the mea- 



^ The performance of a least squares technique can be 
greatly improved by imposing physical limits on the fitted 
parameters. For example I restrict the inclination to lie 
in the range 0° -90° , and the thickness of the gas layer 
to be positive. A standard LSQ routine (e.g. MRQMIN, 
Press et al. 1990) can be easily modified to include phys- 
ical limits in the step where a new guess for the solution 
is calculated: if the new guess exceeds a physical bound- 
ary, then the new guess is placed halfway between the old 
guess and the physical boundary. 



sured apparent widths in the lower velocity chan- 
nels provide enough additional information to 
simultaneous determine the flaring and inclina- 
tion. Furthermore, since many azimuths (chan- 
nels) are used, the effects of local irregularities 
are reduced. 

For less inclined systems, inclination and thick- 
ness can be determined from each channel indi- 
vidually using 



Ad{R; 



chan ) 



COS^(i?;Khan)-2^^(^.^^^^^) 

and 

FWHMl{R-V,^^^) 



(7) 



sin'' i{R) 



[FWHMl^iR; Khan) - FWHM^ + 



[2.35H^c(i?;V'chan)cOSi(i?)]' 



(8) 



where Ad is the projected separation of the two 
"wings" of the "butterfly" . From the inclination 
and thickness values determined from the indi- 
vidual channels the weighted averages are cal- 
culated, where the weights are the inverse of the 
square of the errors. The errors in the right hand 
sides of Eqs. (|^) and (|8|) include the errors on d, 
FWHMobs, Yc, Wc and t. 

Both the intrinsic thickness of the gas layer 
and the planar width of the CSDD increase ap- 
proximately linearly with radius: the former be- 
cause of the 1/R'^ DM-halo density distribu- 
tion (Paper I), the latter because in the wedge- 
approximation of the CSDDs the path through 
this wedge is proportional to R. Simple ana- 
lytic considerations (Oiling &; van Gorkom 19950 
; Oiling 1995b, Chapt. 4) indicate that the ratio 
of these two projected widths depends strongly 
on inclination and only weakly on the proper- 
ties of the galaxy. For a typical channel, the 
projected gas layer width exceeds the projected 
CSDD width for inclinations larger than about 
60 degrees. This inclination cutoff is confirmed 



^ Notice a typo in their Eq. (4): (Jz{R) has to be replaced 
by a^{R) X R 
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using "observations" of simulated galaxies seen 
at various inclinations (Fig. 11; Oiling 1995b, 



Chapt. 4). Thus it is expected that for galaxies 
with inclinations exceeding 60° the flaring, 

and hence the DM-halo flattening, can be deter- 
mined. 

5.2. The Thickness of NGC 4244's Gas 
Layer 

For all major axis positions (x) the emission 
profile perpendicular to the major axis are inves- 
tigated. Gaussian fits to these profiles yield the 
peak and width (Fig. 12) for all positions. All 



the thickness measurements presented in this fig- 
ure, as as well as the odd channels in between, 
are used in the determination of inclination and 
thickness of the gas layer. The peak of the fitted 
Gaussians at each radius is plotted in the two 
outer panels, the apparent thickness in the in- 
ner panels. A fit is only performed if there are 
at least three independent points with a signal 
to noise ratio > 4 in the profile. The SW-NE 
asymmetry in the surface density profile (Fig. |6|) 
at 5.5 kpc and 6-8 kpc can be seen in many of 
the individual channels, indicating that the SW- 
NE asymmetry of the radial surface density dis- 
tribution as determined using the Abel-inversion 
technique extends over a large range in azimuth. 

Because inclination and thickness are calcu- 
lated rather differently in the nearly edge-on case 
[Eq. (^)] and the not edge-on case [compare 
Eqs. (0) and (^)], stringent criteria are ap- 
plied to avoid spurious two-component fits. A 
reliable two-component fit must satisfy the fol- 
lowing criteria: a) the amplitudes of the com- 
ponents may not differ by more than 60%, b) 
the components must have equal widths, c) the 
separation between the components must exceed 
M AX {FWH Mb, 0.77 x FWHMcmp), and d) its 
goodness of fit must exceed the goodness of fit 
of the one-component fit. Generally, the pro- 
files are well represented by a single Gaussian, 
so that, the thickness and inclination are deter- 
mined using Eq. (|6|). The results are presented. 



for the two sides of the galaxy, in Fig. The 
radial variation of thickness in the upper panel, 
the inclination in the lower panel. The plot on 
the left-hand side results if the inclination is a 
free parameter, while fixing the inclination leads 
to the thickness determinations plotted on the 
right. 

Using the observational errors only, the re- 
duced values of the fits [to Eq. (^)] are of the 
order of 3.5 (it 1.5) for radii smaller than ~ 8 kpc 
and increase sharply to ~ 13 (it 7) beyond, indi- 
cating that either the model is bad or that the er- 
rors have been underestimated. In Appendix |G| I 
show that non-circular motions contribute to the 
error budget of Eq. (P) as well. The inclusion of 
errors associated with non-circular motions (~ 2 
km s^"*^ for R < 8.5 kpc and ~ 8 km s^"*^ beyond) 
decreases the reduced to values around one. 
The onset of the warp around 8.5 kpc suggests 
that modest non-circular motions are associated 
with the warp. Other explanations are possible 
however: for example the surface density distri- 
bution may be strongly asymmetric in azimuth 
so that the H I emission arises from either near 
or far side rather than from both sides, or the 
warp may be asymmetric like the warp of NGC 
4565 (Rupen 1991b). 

The beamsize corrected widths inferred from 
the high (lO'fl) and intermediate (38" ) resolu- 
tion data sets are generally equal to within the 
errors, which indicates that the width measure- 
ments are not hampered by resolution effects. It 
is likely that the non-circular motions inferred in 



the inner 5 kpc (Sec. 4.2., Fig. |9|), which might 
result from the "spiral structure" as seen in the 
H I surface density distribution (Fig. ^), are not 
limited to the regions close to the major axis. 
Because the determination of the thickness of the 
gas layer is less accurate in the presence of non- 
circular motions]^, it is not advisable to use this 



^'^In fact, in Appendix |c] I show [Eq. (C21)] that the errors 
in the derived thickness are proportional to the magnitude 
of the non-circular motions and to the inverse of the ro- 
tation speed, which is rising in the inner 5 kpc. 
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region to draw inferences about the mass-to-light 
ratio of the stellar disk. The derivation of the 
halo's flattening, which uses the flaring data be- 
yond 8 kpc, is not affected. 

Another problem area may be the 7 - 8.5 
kpc hump on the south-western side, where the 
thickness of the gas layer is 70% larger, and the 
gaseous surface density distribution (Fig. |^ 50% 
lower than on the north-eastern side. This re- 
gion is adjacent to a region where the stellar sur- 
face density is possibly depressed a factor of ten 
(Sec. [2] and Fig. |2|). In the outer part of the 
optical disk the self-gravity of the gas can con- 
tribute significantly to the vertical force (see Pa- 
per I, Fig. 10). In the case of NGC 4244, models 
which include the gaseous self-gravity have gas 
layer widths 50% smaller than models which do 
not, between 6 and 10 kpc. Therefore I plot the 
north-eastern rather than the average values in 
the final flaring curve (Fig. 14) from from 6 to 
8.5 kpc. 

From these flaring measurements taken at dif- 
ferent resolutions and sensitivities I extract two 
flaring curves, one for the varying inclination 
case, the other for the fixed inclination case (Fig. 
P^). For radii between 6 and 10 kpc I use the flar- 
ing data extracted from a 357-pc resolution data 
set. The high (177 pc) and intermediate (665 pc) 
resolution flaring results are used for the regions 
smaller than six and larger than ten kpc, respect- 
ively. The "errors" are either the weighted error 
on the average of the two sides, or if larger, half 
the difference between the two values, so that the 
systematic differences are reflected by the extent 
of the error bars. 

5.3. Inclination effects 

Because the inferred short-to-long axis ratio of 
the halo is proportional to the square of the true 
width of the gas layer (Paper I, Paper III), the 
thickness of the gas layer has to be measured re- 
liably. It must be stressed that the inclination is 
a critical parameter for determining the intrinsic 
thickness of the gas layer. This is so because, at 



inclinations around 85° a large part of the appar- 
ent width arises due to the projected separation 
of the CSDD wings (Sec. 5.1.), which changes by 



40% between 85° and 83° . At these inclinations 
the DM-halo shape can be determined reliably 
only if the inclination is accurately known (Pa- 
per III). 

5.4. Non-circular Motions 

The assumption of circular orbits is basic to 
the thickness determination technique presented 
in this paper. Closed orbits of ellipticity e^^p ex- 
ist in elliptic logarithmic potentials, for which the 
radial and tangential velocity can be found ana- 
lytically [e.g., Rix & Zaritsky 1995, their Eqs. 
(7) and (8)]. As a result, the difference be- 
tween the observable recession velocities of an 
elliptical disk and a circular disk have a period 
of 7r/3 and an amplitude of et^^pFrot- With 6^ 
the direction of the major axis of the poten- 
tial, the observable velocity difference between 
two points at {R, 6) and {R, —9) equals zero and 
(2e<^^pFrot sin (26*) sin(6' -|- 20^)) for a round disk 
and an elliptical disk, respectively. For elliptici- 
ties as reported by Rix & Zaritsky (0.045^002)' 
the effects of elliptical orbits are smaller than 9 
kms~^, so that the shape of the CSDD maps, 
and hence the determination of the flaring is not 
affected much. The recession velocity at the ma- 
jor axis equals Vgys ± ^rot(^) {1 + e</,,p cos {20^)]. 
Given the predicted and measured upper limits 
to e^,p, (< 0.1, Dubinski 1994; and < 0.075, Rix 
& Zaritsky 1995), the measured 15% decrease in 
rotation speed cannot be explained by a rapid 
radial variation in ellipticity or location of major 
axis of the dark matter distribution. 

6. DISCUSSION 

The thickness of the gas layer has been mea- 
sured for only a few galaxies: The Milky Way 
(e.g., Kulkarni et al. 1982; Knapp 1987; Wouter- 
loot et al. 1990; Diplas & Savage 1991; Mer- 
rifield 1992; Malhotra 1994, 1995), NGC 891 
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(Sancisi & Allen 1979; van der Kruit 1981; Ru- 
pen 1981a), NGC 4565 (Rupen 1991a), and M31 
(Braun 1991). The rapid variations (up and 
down by a factor of 2 - 4 within 500 pc) in the 
thickness of the gas layer of NGC 891 (Rupen 
1991a) and the 7-kpc hump found in the flaring 
curve of NGC 4244 (Fig. H) are probably caused 
by spiral arm streaming motions which change 
the shape of the emitting regions of the channels 
(see also Appendix |^). Inside the stellar disks, 
the thickness (FW H M) of the gas layer increases 
gradually from ~150 pc to ~800 pc at the edge 
of the stellar disk, similar to the flaring of NGC 
4244 as determined in this paper. 

Inside the optical disk {R < 10 kpc), the radial 
variation of the dispersion is quite similar to that 
observed in other, more face-on, systems (van 
der Kruit k Shostak 1982; Dickey et al. 1990; 
see Kamphuis 1993, Chap. 12, for a review). Be- 
yond the optical disk {R > 10 kpc), Ugas increases 
somewhat (from 8.5 to 9.2 km s~^), similar to the 
dispersion increase (from 6 to 8 km s^^) seen in 
NGC 628 (Kamphuis 1993, Chap. 12, Fig. 1). 
In the other two systems where measurements of 
the dispersion beyond the optical disk exist, CTgas 
declines very slowly, to 6 and 7 km s~^ (Dickey 
et al. 1990; Kamphuis 1993). Thus, the radial 
variation of velocity dispersion in the gaseous 
medium of NGC 4244 is not unusual. 

Rather than determining the best fit M/C 
and dark-halo's core radius and central density I 
use the procedure outlined in Paper I to find the 
range in stellar mass-to-light ratios which yield 
acceptable fits to the observed rotation curve 
(the halo's flattening is not a parameter in this 
procedure). I then calculate the reduced val- 
ues on a grid which covers this range in M/C . 
Recall that for given observed rotation speeds 
at 2.3/iR, (14ot(2.3/iR) and 8/ir, (V"rot(8/iR)) and a 
choice for the stellar M/C the dark halo param- 
eters are uniquely determined (Paper I). How- 
ever, because measurement errors on the rotation 
curve, different values values for 14ot(2.3/iR) and 
^ot(8/iR) are acceptable. I adopted values which 



differ by less than l-a from the adopted rotation 
curve which yield the lowest values. In Fig. 
15 I present the so found 'best fit' and the it l-a 
mass models. The mass-to-light ratio is rather 



well determined {M/Cb = 3.71 



+2.6\ 



The up- 



per limit coincides with the maximum-disk value 
(van Albada &: Sancisi 1976), for which the peak 
of the rotation curve due to the stellar disk equals 
the observed peak rotation speed. 

The shape of the vertical distribution of the 
stars determines the inferred mass-to-light ratio 
of the stellar disk (van der Kruit 1988). Using a 
sech^ vertical distribution, the inferred mass-to- 
light ratio is twice larger than when using a verti- 
cally exponential distribution (for edge-on galax- 
ies). As the M/Cb value found here (6.3) for 
the maximum disk case is about twice as large as 
typical for a maximum-disk Scd galaxy (Broeils 
1992, Chapt. 10) the vertical distribution may 
be exponential rather than sech^. Additional in- 
formation (e.g., K-band imaging) regarding the 
vertical structure of the stellar disk is required to 
determine M/C more accurately. 

Beyond the optical disk, the dynamical mass- 
to-light ratio increases to a value of 14.9 Mq/Lq 
at the last measured point (16 kpc), or half that 
value for the vertically exponential distribution. 
Neither value is unusual (Broeils 1992, 1995). 
Only a small amount of dark matter is present: 
1.4 times the luminous mass for the maximum 
disk case, and 3.6 times the luminous mass for 
the half-maximum disk case (independent of the 
vertical distribution). As a result of the limited 
extend of the measured rotation curve, the dark- 
halo parameters are ill constrained: to within a 
factor 50 for the core radius and 900 for the cen- 
tral density. 

7. CONCLUSIONS 

The rotation curve of the almost edge-on ga- 
laxy NGC 4244, determined from high resolution 
H I spectra, rises slowly to about 100 km s~^at 
5 kpc, remains flat till the edge of the optical 
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disk at 10 kpc and decreases in Kepler ian man- 
ner till the last measured point at seven optical 
scale lengths (14 kpc). The mass-to-light ratio of 
the stellar disk lies between 50% and 100% of the 
maximum-disk value, the core radius and central 
density of the dark halo could not be determined. 

I developed a new technique to determine in- 
clination and flaring for large H I rich galaxies 
seen at inclinations larger than about 60° . Kine- 
matic information is used to disentangle flaring 
and warping. The thickness of NGC 4244's gas 
layer increases gradually from ~ 400 pc at 2 opti- 
cal scale lengths to ~1.5 kpc at the last measured 
point. 
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APPENDICES 



A. DETERMINING THE ROTATION 
CURVE AND VELOCITY DISPER- 
SION OF AN EDGE-ON GALAXY 

In a manner similar to the tangent-point method 
employed for the Milky Way (e.g., Burton & Gor- 
don 1978; Celnik et al. 1979), the gaseous veloc- 
ity dispersion and the circular velocity can be 
derived from the extreme velocity tails of spec- 
tra taken on the major- axis of an edge-on ga- 
laxy. In this section I investigate how gradients 
(along the line-of-sight) in surface density, ro- 
tation speed and velocity dispersion affect this 
derivation. Small gradients occur naturally in 
the azimuthally symmetric galaxy models con- 
sidered in this paper, while strong gradients can 
occur near spiral arms. Gradients perpendicu- 
lar to the line-of-sight, beam smearing, are also 
addressed. 

For a galaxy, in which all gas moves at ex- 
actly the circular velocity (K:), gas on the major 
axis moves at the local rotation velocity. Other 
positions have lower projected velocities, Vp [Eq. 
(^]. In reality, the velocity dispersion of the gas 
broadens the sharp cutoff at Vc{R) and moves 
the peak of maximum emission closer to the sys- 
temic velocity. The spectrum, /(x,^, d;5), at a 
point|^ (x, d) on the sky depends on the surface 
density distribution of the gas, the distribution 
perpendicular to the plane, the gaseous veloc- 
ity dispersion, the rotation curve, and the chan- 
nel response function Siy' ]V, /IV), defined by 
Eq. (pi). After summing the data along the mi- 
nor axis, the dependence on minor axis position 
and the vertical distribution disappear so that 
the spectrum (or "XV plot") is given by: 



I{x,V-S) 



K 



dy 



2ttJ " a{y) 
dV S{V';V,^V) e(-i[(^'-^p(?'))/'^(2^)l'),(Al) 

with K a proportionality constant related to the 



^^The coordinate system is defined in Sec. 
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choice of units, y = \/^ 
V^{y) = K(i?)x 



and 



~ K(x) + ^(ax-yc(x)) . (A2) 

Here Vc{R) ~ Vc(x) + a {R — x), and a the 
gradient in the rotation curve (km s~^/kpc). A 
point somewhere along the hne of sight whose 
projected velocity is a km s~^ lower than the ro- 
tation speed will hardly influence parts of the 
spectrum at velocities larger than Vc{x). Thus, 
the spatial integral in Eq. ( |A1| ) can be limited 
to a maximum, iymax so that only gas "close" 
to the major axis will contribute significantly to 
the extreme velocity tails (EVTs). Eq. ( [A2| ) pro- 
vides an estimate for yma.x- 



2 a{x) 



ax 



Vc{x) 



(A3) 



For NGC 4244, typical values for e^{= ymax/x) 
decrease from 0.85 at 1 kpc to 0.33 at 12 kpc. 
The velocity dispersion and the surface density 
can be expanded around R = x as well: cr{y) ~ 
a{x) (l + ^ bxe'^^, and S(y) « S(x) + \ cx<? , 
where h is the fractional gradient of the velocity 
dispersion (kpc^^), and c the gradient in surface 
density (M0pc~^/kpc). 

For spectral line data obtained with the VLA 
which are not Hanning smoothed (as is the case 
for the observations presented in this paper), a 
channel of width Al/, centered at V ^ is sensi- 
tive to emission at other velocities V as well. 
The channel response function is given by (Rots 
1990): 



5(y';y,Ay) 



sin(^(y'-y)/Ay) 



(A4) 



7r(y' - V)I^V 

The strong sidelobes of S are reduced by more 
than an order of magnitude when integrated over 
the width of a channel (the 1**, 2"'^, 3'"'^, and 
4*^^ channels have an effective sensitivity of 100%, 



-|-8.6%,-1.9%, and -hO.83%, respectively). Limit- 
ing the velocity range (to zt^AF) of the second 



integral in Eq. ( |Al| ) and taking S to be unity 
over this interval allows for an analytic approxi- 
mation of I(x, V] S), which can then be used to 
gauge the importance of line-of-sight gradients of 
the variables. 

Approximating the Gaussian term in Eq. ( |Al| ) 
by a parabola and integrating the spatial integral 
up to em, the spectrum can be written as: 

Em K Ay 



I{x,V) 



2Txa 



'2i <T ; X 




+ 



c 



(A5) 



where all galaxian properties are evaluated at ra- 
dius R = x, and with W = V — Vc{x). The multi- 
plicative term in Eq. ( |A5| ) is a Gaussian function 
of dispersion a. The first term within the curly 
brackets results from material on the major axis, 
the second term arises due to material close to 
the major axis, while the last three terms orig- 
inate due to radial gradients in rotation speed, 
dispersion and surface density, respectively. Ex- 
cept in the inner 6 kpc, where beam smearing is 
important as well, Eq. ( [A5| ) represents the mea- 
sured EVTs well. A beam smearing correction 
could be obtained by from an approximation of 
the convolution of the spectrum (|A5| ) with the 
beam. However, as this leads to quite ugly equa- 
tions, another approach is followed. As we have 
seen above, the EVTs of an XV plot contain in- 
formation about the intrinsic rotation curve and 
velocity dispersion. A Gaussian fit to such an 
EVT yield a value for the position and width of 
this GaussianP|. Line-of-sight and beam smear- 
ing effects cause these values to deviate from the 



^■^In practice multiple components are fitted to the full spec- 
trum such that the EVT is dominated by one component. 
Because of the additional components (more degrees of 
freedom) the errors are overestimated. The AIPS tasks 
SLICE and SLFIT were used to perform the Gaussian 
fits. However, the version of these tasks which are part of 
the standard AIPS release re-sample every slice to a 512 
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intrinsic rotation speed and velocity dispersion. 
The magnitude of these deviations can be em- 
pirically found by comparing the values obtained 
from fits to a simulatedP] XV plot with the pa- 
rameters used in the calculation of the simulated 
XV plot. I find: 



C^XVilO 



gas 



(0.9-Fexp[-i?/4]) 



(A6) 



,(A7) 



^ (0.9 exp[-2.5] + exp[(i? - 10)/17]) 

for (R < lOkpc) and {R > lOkpc) respectively. 
The subscripts "10", "38", and "55" refer to the 
beamsize of the data set for which the correction 
is appropriate. The "XV" subscript refers to the 
value as fitted to the extreme velocity envelope of 
the XV plot. Note that the true gaseous velocity 
dispersion is always smaller than the dispersion 
of the extreme velocity envelope. As these correc- 
tions depend on the input parameters (e.g., beam 
smearing is insignificant for a flat rotation curve) , 
the input parameters for the calculation of the 
simulated XV plot should be close to the intrin- 
sic rotation and dispersion curves: this problem 
is solved by iterating several times. An estimate 
of (Txv can also be obtained by counting the num- 
ber (A^o,p) of velocity channels between the zero 
and peak emission level of the EVT. This esti- 
mate is calibrated by generating Gaussian dis- 
tributions with a range of dispersions centered 
at various locations within a velocity channel: 
(axv = (2.4iVo,P - 0.6 ± 1) x kms"^). 
Since A^o,P as determined from the data is gener- 
ally equal to the value determined from the best- 
fit simulated XV plot, I conclude that the ve- 
locity dispersion (and rotation curve, see below) 
as determined from the Gaussian-fit method is 
reliable. 

Because these relations depend on beam size, 
channel width, surface density distribution, and 



grid, as a result of which the calculated errors are mean- 
ingless. I therefore wrote special versions of SLICE and 
SLFIT which do not re-sample the data. 
Using the procedure outlined in Appendix p]. 



the rotation curve, these relations are specific to 
the observations of NGC 4244 presented in this 
paper. The errors, 3a, about 20% and 5% in the 
inner and outer part of the galaxy, are largest in 
those regions with strong surface density gradi- 
ents. 

The corrections required to determine the ro- 
tation curve from the XV-plots are given by: 

VrotAO (^) ~ ^XV;10 + 

(1.3 X exp[-i?/4.0] + 0.35) x o-xv;io (A8) 

Kot;38,55(-R) ~ VxV';38,55 + 

(1.0 X exp[-ii/4.0] + 0.35) X f7xV;38,55. (A9) 

In regions where the surface density changes 
rapidly, the smearing correction is least reliable: 
differences between input and output rotation 
curves of order 3-5 km s~^ occur in the inner 5 
kpc of the simulated XV plots. 

These correction are then applied to the fit- 
ted rotation and dispersion values obtained from 
the data. The resulting dispersion and rotation 
curves are shown in Figs. ^ and ^, respectively. 
At some locations, the data EVTs are very dif- 
ferent from the simulated EVTs, while adjacent 
positions have similar EVTs: these locations are 
likely to be affected by local irregularities and 
have been discarded. 

This procedure does not work perfectly, not 
even for the azimuthally symmetric model ga- 
laxy used here. Both the XV-derived rotation 
curve and the XV-derived dispersion curves (af- 
ter applying the corrections specified above) de- 
viate systematically from the input curves. Be- 
tween 1 and 2 kpc and between 4.5 and 5 kpc 
the XV-derived rotation curve falls ~8 km s~^ 
(~12%) and ~ 3 km s~^ (~3%) below the in- 
put curve. These negative differences can be 
clearly attributed to the strong density enhance- 
ments at slightly larger radii, which are thus 
seen at lower projected velocities: on the "in- 
side" of strong H I density enhancements the 
XV-velocities are decreased and the true rota- 
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tion speeds are underestimated. Too large a ro- 
tation speed (3 km s~^) is derived in the region 
where the H I distribution strongly peaks (Fig. 

between 3.5 and 4.5 kpc. This H I peak be- 
haves like an 5 function on the major axis, in 
which case the XV-velocity directly corresponds 
to the rotation speed: applying the dispersion 



correction [Eq. (A8)] thus leads to an overesti- 
mate of the rotation speed. In determining the 
rotation curve of NGC 4244 (§|42l) I have not 
tried to compensate for the surface density effects 
in detail since these H I density peaks may well 
be spiral arms, in which case the derived rota- 
tion speeds may include some spiral arm stream- 
ing motions as well. Instead I have followed the 
upper envelope of both sides and in the regions 
where the model calculations indicate that too 
low a rotation speed is inferred from the XV- 
plots I increased the adopted rotation curve such 
that is a smooth junction of the adjacent regions. 
Here the errors are accordingly increased. 

Because beam smearing effects are not impor- 
tant beyond 6 kpc, Eq. ( |A5D can be used to 
investigate the dependence of the shape of the 
EVTs on gradients in Vc{R), T.{R), and a{R). 
Since for a given range in y (or em) only a small 
range in R is sampled, gradients in the rotation 
curve are unimportant. Gradients along the line 
of sight in surface density and dispersion cause 
changes in Vxv and fixv- Using Eqn. ( |A5|) and 
determining Vxv and axY as a function of b and 
c, I find the following offsets (in km s~^): 



AFxv « (-3 ± 0.7) b , 



(AlO) 



AcTxv ~ (1-8 ± 0.8)6 + (0.84 ± 0.12) ( ^ ) (All) 



for (c < 0), while the c-dependence in ( All| ) is 
absent for c > 0. Where the A-symbol repre- 
sents the difference between the "no gradients" 
and the "strong gradients" XV-parameters. Be- 
yond the optical disk, typical values for c/E and b 
are -0.5 and 0.05, respectively. The smearing cor- 
rections [Eqs. (A6) through ( |A9| )] take these gra- 
dients into account. From Eq. ( |A11 ) it follows 



that in order to over estimate axv (and hence 
cr) large positive line of sight gradients in the dis- 
persion are required. With axv of the order of 
10 km s"\ a 10% over estimate is possible if the 
dispersion along the line of sight increases to 15 
km s~^ over a distance of 1 kpc. If such a strong 
(T-gradient would exist in the radial direction as 
well, it would certainly be noticed. Very large ve- 
locity dispersion gradients have to exist to affect 
Vxv significantly. I conclude that the rotation 
and velocity dispersion curves can be determined 
reliably from XV-plots of edge-on galaxies. 

B. SIMULATING VLA SPECTRAL 
LINE OBSERVATIONS 

In order to facilitate the interpretation of 
H I spectral line data I developed software which 
simulates VLA H I observations of galaxies. The 
accuracy of the determination of rotation and 
velocity dispersion curves from XV plots, the 
beam-smearing corrections, the thickness-and- 
inclination determination technique, can all be 
gauged by comparing the "observations" of such 
a simulated galaxy with the input parameters. 
Input parameters to this program (VLASIM) 
are: the observing geometry, S(i?), VrotiR), and 
CTgas • The vertical density distribution of the gas 
layer is calculated employing a physical model for 
the potential of the galaxy (the multi-component 
method. Paper I) which yields the gaseous vol- 
ume density. The rotation speed is assumedf^ 
to be constant with height above the plane. 
VLASIM then integrates the amount of hydro- 
gen along lines of sight from the observer through 
the galaxy, and distributes this amount over the 
appropriate channels [using the channel sensitiv- 
ity function Eq. ([A4|)]. When the line of sight 
integrations are done, the spectral line cube is 
smoothed to the requested resolution followed by 
the addition of realistic noise (if so desired). 



^''For several disk-halo models, rotation curves at different 
z-heights were calculated, all showing that the rotation 
speed does not vary significantly (< 5 %) in the region 
where gas is found (| z | < 3FW H M^.gas)- 
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C. THE EFFECTS OF KINEMATICAL 
DISTORTIONS 

As discussed in Sec. knowledge of the ro- 
tation curve allows us to determine the intrinsic 
thickness of the gas layer. However, the tech- 
nique is sensitive to deviations from the assumed 
rotation curve. Spiral arm streaming motions 
of several tens of percent of the rotation speed 
(several tens km s~^) are not uncommon in spi- 
ral galaxies (e.g., Visser 1980; Vogel et al. 1988; 
Boulanger & Viallefond, 1992; Tilanus & Allen, 
1993), so that locally the velocity deviates from 
the circular velocity. When a line of sight is along 
a spiral arm, there will be a net shift in the radial 
velocity (type 1), when the line of sight crosses 
a spiral arm, the radial velocity will not be af- 
fected but the dispersion will be larger]^ (type 
2). In this section I determine which parts of 
the galaxy are most susceptible to non-circular 
motions and present a weighting scheme which 
minimizes the effects of non-circular motions: a 
summary is presented in the next section ( |C.l. ). 

The galactocentric coordinates [{x, y, z) or 
(i?, 0, z)] are related to the celestial coordinates 
[(x,d) centered on the galaxy] in the following 
way: 



X = R cos 6 , 

d = R cos i sin -\- z sin i . 



(C12) 
(C13) 



where the position angle of the major axis was 
taken east-west (90° ) and where z is the height 
above the plane. If gas located at (i?, 6, 0) has 
an excess line of sight velocity of 



AV 



Sim 



[T4ot(i?o) cos ^0 - V,ot{R) cos (^0 + 



with respect to a position {Rq, Oq,0), this position 
will be "seen" in the same channel as {Ro,6q,0) 
is. Similarly, a channel Ay km s^^ wide "sees" 
gas in an area delimited by 9q and 9. Expand- 
ing cos (6*0 -|- A^) to second order around 6o, and 



assuming that rotation curve around Rq is con- 
stant, I approximate the above equation and 
solve for AO. 



AV CSC i 



VrotiR 



and 



sin 00 Ae + -cos9o{ABf 



tan 9t] 




2Aycsci cos 01 



'0 



, (C14) 



(C15) 



V,ot{Ro) sm'9o 

The requirement that the two points [(i?0)^0)0) 
and (i?, 9, 0)] have the same major axis distance 
leads [via Eq. ( |C12| )] to the following relation for 
AR {=R-Ro): 



AR 



A9 



1- A9 tan 9, 



tan 6*0 (C16) 



In the lower panel of Figure 16 I present the 
relations, both practically independent of incli- 
nation, for A9 and AR graphically for the case 
where AV = 10 km s"^ and yrot(-Ro) = 100 
km s~^. Although this galaxy model is extremely 
simple, the opening angles as determined for a 
more realistic model behave similarly: the high 
velocity channels (small 9o) are wider than the 
low velocity channels. For 9o < 50° the frac- 
tional range in galactocentric radius is small so 



I thank Ron Allen for pointing this out. 



that the approximate relation for AV ( C14 ) is 
reasonable, whether or not the rotation curve is 
constant around -Ro- 
lf the positions {Rq, 9q, zq) and {R, 0, 0) project 
on the same position on sky, and if due to an 
excess motion they also have the same velocity, 
then their positions in the observer's spectral-line 
cube will be the same. The equivalent vertical 
height of this streaming motion (or, equivalently, 
of a channel AV km s^^ wide) equals zq. Using 
( CIS ), ( C16| ) and expanding sin(0) around 9o I 
find: 

« Ro ( . ) . (C17) 

Vcos 9ota.ni J 

I plotted the projection of zq (= zosini) in the 
second to bottom panel of Figure |l^. The curves 
are labeled with the corresponding galactocentric 
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radii. For i = 85° its value is small compared to 
the projected vertical thickness of the gas layer, 
but its importance relative to the projected ver- 
tical thickness of the gas layer grows rapidly with 
decreasing inclination (oc tan~^ i). 

The systematic errors introduced by stream- 
ing motions are more easily understood when 
considering the corresponding planar width of 
the CSDD, AY{AV) ^ Rq (^), which can 
be approximated as: 



AY{AV) 

Rq 



COS 6*0 (1.81V sin^-^(9o + 



v^Vcos 6*0 



-1 



(C18) 



with V = [(VrotiRo) smi)/{2AV)]. In the region 
close to the major axis, kinematical distortions 
have an effect proportional to the square root 
of the fractional velocity deviation (2"'^ term), 
in the lower velocity channels the proportion- 
ality is approximately linear (1** term). Thus, 
the channels close to the major and minor axes 
have long pathlengths through the galactic disk 
and are hence easily affected by kinematical or 
surface brightness irregularities (Fig. 16, third 
panel from below). The edge-channel itself is 
even wider as the two sides (at =b^o) merge to- 
gether. 

Streaming motions of type 1 cause the centers 
of the CSDD to shift along the line of sight, typi- 
cally by ±AY [see also Eq. ( ^201 )]. In the almost 
edge-on situation, the two sides of the CSDD, as 
projected on the sky, are overlapping so that the 
projected planar width [Eq. (^)] and its error 
are dominated by the Yc term. Visser (1980) 
presents model velocity fields for M 81 which in- 
clude spiral arm streaming motions (his Fig. 10). 
The equivalent CSDD map results after deprojec- 
tion of the velocity field (see also Sec. 5.1.). From 



these velocity fields it is clear that the centers of 
the CSDDs shift substantially due to spiral arm 
streaming motions: outwards on the outside of 
the arms and inwards on the inside of the arms. 
The strong kinks follow the spiral arm and would 



produce sudden changes in the separation of the 
near and far sides of the CSDDs and in the width 
thereof, over a large range in azimuth. If stream- 
ing motions are present, but the resulting CSDD- 
kinks are not taken into account in the analysis, 
the inferred widths may show strong gradients 
when crossing a spiral arm. Such may be the 
reason for the large thickness gradients seen in 
NGC 891 (Rupen 1991a) and NGC 4244 (the 7- 
kpc hump. Sec. |5.2.D . 

At lower inclinations, when the two sides of 
the CSDD are separated, streaming motions of 
type 2 affect the width determination [Eq. ( P21| )], 
while the errors on Y (type 1 streaming motions) 
hamper the determination [via Eq. (|^)] of the in- 
clination. This is illustrated in the top panel of 



Figure 16, where I plot the inclination error re- 
sulting from a velocity shift of 10 km s~^: 



Ai 



cos I 



sin^ ^ 



A6 



(C19) 

sm fo cos c^o 

Thus, small streaming motions occurring close to 
either the minor or the major axis can cause large 
errors in the inferred inclination. Furthermore, 
the inclination errors are strongly inclination de- 
pendent: Ai{i = 65° ) 6Ai{i = 85° ). 

C.l. Summary 

I will use the above determined systematic er- 
rors estimates in the thickness & inclination de- 
termination process (Sec. ^.Ij ). The effect of spi- 
ral arm streaming motions is to either shift the 
center of the CSDD along the line of sight or 
to increase the apparent velocity dispersion, the 
former having more severe effects. To our advan- 
tage, the expected systematic errors scale with 
the width of a CSDD wing (calculated for a given 
gaseous velocity dispersion figas), 



6Yc ' 
and 



AY{AV) 

Ay(cJgas) 



X Fiy™c(cTgas) ,(C20) 



FWHMc{(J„ 



(AY{AV)_\ 
Y^[AY{a,.,)) 



1 + 



(C21) 
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are easily calculated from the CSDD model maps 
and Eq. ( |C18 ). These systematic errors increase 
towards lower inclinations and are proportional 
to (sini sin^cos^)"^: the regions around both 
principal axes are most easily affected by stream- 
ing motions. 

D. OPTICAL DEPTH EFFECTS 



In the Milky Way the optical depth in the 
H I line is substantial within a few degrees from 
the galactic plane (e.g., Dickey et al. 1983). Typ- 
ically, the total emission along the line of sight 
is under estimated by 15%. Similar values are 
found for edge-on external systems (Haynes k. 
Giovanelli 1984). We thus expect similar values 
for NGC 4244. Braun & Walterbos (1992) esti- 
mate that the cold interstellar medium, which is 
mostly responsible for the absorption (see below), 
has a volume filling factor of about 1 - 8%, for 
the Milky Way and M31, respectively. This cold 
material forms a layer two to three times thinner 
than the warm H I (e.g., Kulkarni & Heiles 1988; 
Braun & Walterbos 1992). For NGC 4244, we 
thus expect that emission profiles perpendicular 
to the major axis might be affected by absorp- 
tion on the midplane, resulting in larger apparent 
widths. However, from the peak surface bright- 
ness and apparent width measurements (Fig 



12) 



such effects are not obvious: profiles with a high 
peak column density do not seem to be broader. 
A more detailed, albeit idealized, calculation of 
the optical depth is given below. 

The optical depth in the H I line for a uniform 
medium of length / and number density n is given 
by Braun & Walterbos. Their Eq. (7) reads: 



170- 



n I 



ts,100 



Ay 



(D22) 



with Tg^ioo the spin temperature in units of 100 
K, and the width of the channel in km s^^. 
In Sec. ^ the surface density, the rotation curve 
and the velocity dispersion were derived assum- 
ing low optical depth. Those derived quanti- 
ties can be used to evaluate the quantities in 



Eq. ( p22 ) and check the validity of the as- 
sumption of low optical depth. The CSDD maps 
provide the H I surface density for which the 
projected velocity equals the observing veloc- 
ity (ib^Ay). Investigating these CSDD chan- 
nel maps (Fig. IC) we find typical peak col- 
umn densities of 0.5 - 1 M0pc~^ per channel 
for galactocentric radii smaller than 7 kpc, and 
much less beyond. In edge-on geometry, the sur- 
face density distribution along the line of sight 
is a double Gaussian [Eq. (§)] with a typical 
width {FWHMp) for each component of 3 it 1 
kpc. The H I volume number density is given 
by: n = 5.5 IQ'^ {T.i/FWHM^^i) cm'^, where 
the units for surface density and thickness are 
M0pc~^ and kpc. The optical depth can then 
be rewritten as: 

Si FWHMr, 3 
T = 0.156—^—— — , p23) 

where a channel width of 6.2 km s~^ was used. 
Typical values for the spin temperature of the 
cold component of the interstellar medium range 
from 100 - 175 K (for the Milky Way and M31, 
respectively; Braun & Walterbos 1992). Smaller 
values are found in individual clouds (Spitzer 
1978, p. 48, and references therein). Taking the 
model widths as representative for the thickness 
of the gas layer (~ 250 pc. Fig. 14), and assum- 



ing that the cold medium fills the entire volume 
we find r ~ 0.33 (30% attenuation). If the cold 
medium is clumpy with a low filling factor, as 
in the Milky Way, we might expect 3 clouds per 
kpc with an optical depth > 1 (Kulkarni & Heiles 
1988). In an exactly edge-on geometry this den- 
sity of clouds would correspond to a large opti- 
cal depth. However, at an inclination of 85° the 
path through the galaxy (~ 3 kpc) corresponds 
to a change in z of ~ 260 pc (about three times 
the thickness of the cold gas layer), so that the 
number of absorbing clouds per line of sight is 
reduced. Because the pathlengths FWHMp are 
comparable for all velocity channels, the optical 
depth effects are comparable in all channels (in 
the inner 7 kpc). 
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Beyond 7 kpc the surface densities are much 
smaUer so that optical depth effects are not im- 
portant in the derivation of the shape of the DM 
halo (Paper III). 
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Fig. 1. — This image of NGC 4244, taken on a IllaJ-plate, was kindly made available by Piet van der 
Kruit. See van der Kruit & Searle (1981a) for details. The weak dustlane, barely visible on a TV monitor, 
could not be reproduced in this greyscale image. The contours ( [ 2 4 8 16 ] x 6 lO^'' cm~^) represent 
the total H I column densities at a resolution of 10"xlO" (175x175 pc). A greyscale representation of the 
total H I map is presented in Fig. |3[ 

Fig. 2. — In this figure I present four radial light profiles, each again composed of three profiles: SW 
(dotted line), NE (dashed line), and average (full line). For clarity, the profiles have been offset by a 
constant C. From top to bottom the profiles are derived from: a) all light (C=0.69), b) all light within 
420 pc from the major axis (C=-2.81), c) all light more than 420 pc "above" (east of) the midplane 
(C=-6.31), and d) all light more than 420 pc "below" (west of) the midplane (C=-9.81). Extinction is 
not important in the determination of the radial scale length (/ir) because all profiles are consistent with 
the same radially exponential distribution, with /ir=2.0 kpc (the filled circles). 

Fig. 3. — Two representations of the total H I distribution are shown. A low resolution (38"x38" , 
665 pc) contour plot is superimposed on a high resolution (10"xlO" , 175 pc) grey scale image. The 
image has been rotated counterclockwise by 48° (NE is upp, SW is down). The units for the contours 
(0.1 0.4 0.8 1.6 3.2 6.4 12.8) and the greyscales are lO^'^ cm~^. A tickmark along the vertical axis, I' , 
corresponds to roughly 1 kpc. The slight warp sets in around 9 arcmin, the edge of the optical disk. 
Clear features in the high resolution image are the the high surface brightness peaks seen on both sides 
located approximately 6 arcmin from the center. Notice an approximately 2 kpc diameter region of low 
column densities, extending from about +3 to +5 kpc. The edges of this feature seem to extend to larger 
height above the plane. This feature is also clearly visible in the channel maps (Fig. |5|). 

Fig. 4. — The Global Profile, or the total H I flux as a function of velocity, determined from the 
55"x55" resolution data set. The systemic velocity is the mean of the 10%, 20% and 50% points, 244.4 
± 0.4 kms~^. A small asymmetry between the approaching and receding sides can be seen. Typical 
uncertainties in the flux equal 0.025 Jy/channel. 

Fig. 5. — A contour representation of the channel maps at 38"x38" (665 pc) resolution superposed 
on a greyscale image of the channel maps at 28'f3xl0"l, (495x177 pc) resolution. The greyscale levels 
correspond to column densities of [1 1.5 3 6 9] lO'^*' cm^^ (the lowest level is ~ 2.6cr), the line contours 
to [-0.07 0.07 0.1 0.2 0.4 0.8 ] 10^° cm"^ (the lowest contour is ~ 2.5cr). The channels of the approaching 
and receding sides are combined into one frame. The value printed in the upper left corner indicates 
the velocity offset from the systemic velocity (i.e., | Khan — Kys I / AVchan)- In order to visualize the 
vertical structure better, the horizontal scale is expanded. The orientation is as in Fig. the distances 
are relative to the pointing center listed in Tab. ^, 

Fig. |5|b. — The channels which show emission on the receding as well as on the approaching side are 
shown individually here. The same coding of greyscales and line contours is used as in Fig. ^ . 
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Fig. 6. — The gaseous surface density calculated using eqn (|2[ see Sec. on a linear (lower panel) 

and logarithmic scale (upper panel). The full line, the dotted line, and the dashed line represent the 
average, the south-western, and the north-eastern surface density profiles. Like in the Milky Way (Knapp 
et al. 1978; Blitz, private communication), no sharp cutoff is observed in the H I distribution. 

Fig. 7. — The major-axis velocity (XV) plot of NGC 4244. The lowest five contours are derived from 
channel maps at low resolution (55" ), the highest five contours from 10"! resolution channel maps. All 
emission within 2.5 kpc from the major axis was summed. The two sides are remarkably similar, both in 
kinematics and extent (to about 15 kpc or seven optical scale lengths). Only between -|-3 and -|-5 arcmin, 
the north-eastern side shows a marked dip in the XV-plot. This area is readily identified with a lack of 
gas on the major axis in the channels 18 through 20 (see Figure ^). Note that this hole does not "punch 
through" the whole disk. The derived rotation curve is indicated by the filled circles. 

Fig. 8. — The beam-smearing corrected [using Eqs. ( |A6D and (|A7|)] tangential velocity dispersion as 
derived from the high (10"lxlO"l), intermediate (38"x38" ) and low (55"x55" ) XV-plots (bottom, middle 
and top panel respectively). The south-western, north-eastern, and adopted average curves are indicated 
by the dotted line, the dotted full line and the filled circles respectively. The measurements are consistent 
with a constant velocity dispersion of roughly 8.5 km/s. The fat horizontal bars in the lower right corners 
indicate the size of the beam. 

Fig. 9. — The beam-smearing and velocity dispersion corrected [using Eqs. ( [A8| ) and ([A9|)] rotation 
curve for the two sides of the galaxy. For an explanation of the symbols, see Fig. ^. 

Fig. 10. — This figure shows a representative set of CSDDs. The channel numbers are shown in the top 
right corner of each panel. As the CSDDs are symmetric with respect to the major axis (the horizontal 
axis), only one quadrant is displayed. The vertical axis of each plot represents the distance along the line 
of sight (LOS axis). The contours represent the effective column density a given channel is sensitive to. 
In this coordinate system, the observer would be located at (0, — oo). We clearly recognize the wedge-like 
structure valid for the simple model galaxy, however, the wedges have become wiggly and curve back 
towards the major axis (as a result of the declining rotation curve). For channels close to the extreme 
velocity channel (-19), the gaseous velocity dispersion broadens the wedges significantly, and fills up the 
gap between the near and far side. 

Fig. 11. — The results of simulated VLA H I spectral line observation of a model galaxy like NGC 4244, 
inclined by 72° . The simulated spectral line cube is processed in the same manner as the observations 
of NGC 4244. The recovered inclinations are shown in the lower panel, while the recovered widths are 
presented in the top panel. The filled circles represent the input model widths. No noise was added to 
the simulated spectral line cube. 
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Fig. 12. — In this figure the results of the Gaussian fits to the channel maps are presented. The middle 
two panels show the measured widths (FWHM, beam size corrected, in units of kpc), the outer two 
panels show the fitted peaks (Mqpc"^, beam size corrected). The results obtained from the {W".l x 28"3, 
177x495 pc) channel maps are represented by the open squares, the fits to the intermediate resolution 
(38" , 665 pc) maps by the full line (with error bars). At few positions two component fits could be 
made, these are indicated by filled symbols (filled triangles for the intermediate resolution data). The 
two left-hand panels are determined from the south-western side, the two right-hand panels from the 
north-eastern side. 



Fig. 13. — Here I present the thickness (upper panel) and inclination (lower panel) derived from the 



fitted widths (at lO'fl resolution) presented in Fig. [l^ using the procedure described in Sec. 5.1. . Both 
sides (south-western part on the left) of the galaxy are shown. The fat bar in the lower left corner of 
the top panel represents the size of the beam, horizontally and vertically for the major axis and minor 
axis, respectively. For the panels on the left hand side I allowed the inclination to vary, the panels on the 
right hand side are calculated with a fixed inclination of 84.5°. The filled circles represent the average 
of both sides. The weighting scheme which takes potential errors due to streaming motions into account 
was used (see Appendix [C]). 



Fig. 13b. — Thickness and inclination as derived from the 38" x 38" resolution channel maps. 
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Fig. 14. — The lower panel shows the measured widths (open and filled triangles with error bars, for the 
fixed and varying inclination case, respectively). The observed flaring curve is the average of the north- 
eastern and south-western sides, except in the 6.5 - 9 kpc region, where the north-eastern flaring data 
are used (see Sec. 4^) The 7 kpc hump is probably caused by spiral arm streaming motions (Appendix 
^). The top panel displays the velocity dispersion measurements (full line with error bars). The dotted 
lines are determined from the south-western side, the dashed lines from the north-eastern side (lO'fl and 
38" , for i? < 9 kpc; 38" for R £ (9, 13] kpc, and 55" beyond). The adopted dispersion curve (full line 
with error bars) is the weighted average of these individual curves. 
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Fig. 15. — The top panel shows disk-halo decompositions of the observed rotation curve (squares with 
error bars) into components due to the gaseous disks (H I , crosses; H2, open circles), a component due 
to the stellar disk (vertically sech^ and radially exponential till i?niax=10 kpc, M/Cb =1; filled circles) 
and a DM component. The individual contributing DM components are not shown. The data points are 
independent out to 13 kpc, beyond the resolution equals 1 kpc. The rotation values at 2.3/ir and 8/ir 
are differ slightly (100 and 87 km s~^) from the adopted rotation curve (Fig. |9|) such as to minimize the 
residuals (lower panel), see text. The three disk-halo decompositions ([M/Cb = 6.33, Rc = 13 kpc, /Jh,o = 
1.43 mMe/pc^] full line, [M/Cb = 3.71, = 0.65, ph,o = 233] dotted line, [M/Cb = 3.27, i?c = 
0.28, /9h,o = 1300] dashed line) represent the best fit {M/Cb =3.71, x^=0.864, 32 degrees of freedom) and 
those models with a reduced value 1.0 larger. For the disk-halo decompositions presented here, the 
/m, /si and /?8 parameters (as defined in Paper I) have the values 0.7094, 0.6033, and 0.861 respectively. 
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Fig. 16. — This figure illustrates the possible effects of streaming motions on the apparent widths (see 
Appendix C). All curves presented scale roughly linearly with the ratio of streaming motion and circular 
speed. For CSDDs centered on I plot the angular width (the lowest panel, crosses) and the normalized 
range in radius (filled circles) which corresponds to a streaming motion of magnitude AV km s"^. The 
corresponding apparent width on the sky is plotted in the one but lowest panel. In the second panel from 
the top, I plot the linear distance, in the plane of the galaxy, that the angle A6 corresponds to. Close to 
the major and minor axes, long pathlengths reflect a strong sensitivity to small streaming motions. The 
inclination errors resulting from such streaming motions are plotted in the top panel. 
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Table 1: Optical parameters 



Distance 




3.6 Mpc 


B-band radial scale length 




2.0 kpc 


B-band vertical scale height exponential 


Ze{ = ^Zo) 


0.21 kpc 


edge of optical disk 


-Rmax.SW 


10.0 kpc 




-RmaXjNE 


8.0 kpc 


face-on central surface brightness 


MO,B 


22.4 ± 0.07 mag arcsec"^ 




Lb(0) 


74.6 ± 5 L© pc-2 


diameter of 25 mag arcsec"^ isophote 


Dh 


11 kpc 


total B-band luminosity 


Ltot,B 


1.8 10^ L© , Mb = -17.7 



Note: No correction for galactic extinction was required (E(B — V)=0; Burstein & Heiles 1984). The subscript "B" refers to 
the value as determined from B-band photometry. The face-on central surface brightness is calculated assuming a sech^{z/zo) 
vertical light distribution. For a distribution which is exponential in the vertical direction, L(0) is twice larger. The absolute 
magnitude, the central surface brightness, and the total luminosity, which are derived from the central surface brightness 
and /iR.B , change accordingly (to 21.6 mag arcsec"^ , 149 Lq pc~^ and Mb = -18.4, respectively). Note that Dl^-values 
derived from the sech^ and exponential luminosity density distributions (9.6 and 12.6 kpc) bracket the statistically corrected 
RC2 value (de Vaucouleurs et al. 1976). 



Table 2: Instrumental parameters for the VLA observations of NGC 4244 



Observing dates 


2 May 1989, B-array, 6 hours on source 




30 Nov 1989, D-array, 15 min on source 




29 Nov 1990, C-array, 45 min on source 




21 Aug 1992, D-array, 7 hours on source 


Pointing center 


12^ 14"" 59?9 




38° 5' &'.0 


Velocity of band center 


245 km s~^, heliocentric, optical definition 


Number of channels 


63 


channel separation 


o.z Km s 


Velocity resolution 


6.2 kms"^ 


map making characteristics 


beam sizes: along and _L to the major axis 


lO'.'O X 10".0 


B -|- C -|- D-array, uniformly weighted, tapered 


28'.'3 X 10"1 


B -|- C -|- D-array, uniformly weighted, tapered 


SB'.'O X 38"0 


B -I- C -I- D-array, naturally weighted 


55'.'0 X 55"0 


1992 D-array, naturally weighted 


rms noise in channel maps 




lO'.'O X 10"0 


2.3 mJy beam^i 131 10^^ cm^^ 


28'.'3 X 10"! 


1.9 mJy beam^i 39 10^^ cm^^ 


38'.'0 X 38';0 


0.73 mJy beam"! 2.8 10^^ cm'^ 


55'.'0 X 55"0 


0.76 mJy beam"! 1.4 10^® cm'^ 



Table 3: Parameters derived from the I- 



aadopted(1950) 


12^ 14'" 59!59 


(5adoptod(1950) 


38° 4' 59"7 


Position angle 


-48° 


Inclination 


- 84°5 


I 


11.0 kpc 


Vsys, heliocentric 


244.4 ± 0.3 km s"^ 




(1.34 ± 0.07) 10^ Mq 


Mtot.H I /Lb 


0.74 MqLq-i 


W20 


180.4 kms"^ 


W50 


180.6 kms"^ 



I distribution 



Notes: The adopted center of mass is derived from the total H I distribution (Sec. The H I -radius, j , is the radius 
at which the H I surface density equals 1 Mq pc~^. W50 and W20 are the corrected widths of the global profile at 50% and 
20%, respectively, calculated according to Botinelli et al. (1990). . 



